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Abstract. We employ the filter-design perspective and derive the filter functions 
according to nested Uhrig dynamical decoupling (NUDD) and Symmetric dynamical 
decoupling (SDD) in the pure-dephasing spin-boson model with N qubits. The 
performances of NUDD and SDD are discussed in detail for a two-qubit system. The 
analysis shows that (i) SDD outperforms NUDD for the bath with a soft cutoff while 
NUDD approaches SDD as the cutoff becomes harder; (ii) if the qubits are coupled to 
a common reservoir, SDD helps to protect the decoherence-free subspace while NUDD 
destroys it; (iii) when the imperfect control pulses with finite width are considered, 
NUDD is affected in both the high-fidelity regime and coherence time regime while 
SDD is affected in the coherence time regime only. 
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1. Introduction 

Decoherence of quantum system is inevitable, since coupling of a system to its 
surrounding environment exists all the time. Suppressing decoherence is significantly 
important for a variety of fascinating quantum-information tasks. In particular, 
protecting generated entangled states from decoherence is an important issue that has 
been addressed in different frameworks P-[5].0ne of the most effective techniques to 
combat with the decoherence process is dynamical decoupling (DD) [6l[7], which evolves 
from the Hahn echo [8] and develops for refocusing techniques in nuclear magnetic 
resonance (NMR) P-FTT]. The central idea of DD is first explicitly introduced to 
protect qubit coherence in f6], and it is soon incorporated within a general dynamical 
symmetrization framework |7j. DD schemes operate by subjecting the system of interest 
to suitable sequences of external control operations, with the purpose of removing or 
modifying unwanted contributions to the underlying Hamiltonian. Several DD pulse 
sequences, such as Carr-Purcell-Meiboom-Gill (CPMG) p^llTS] . concatenated dynamical 
decoupling (CDD) [HHIG] and Uhrig dynamical decoupling (UDD) [T7j, have been 
proposed by now. The above DD sequences focus on single-qubit systems, but some 
studies have been extended to multi-qubit systems recently [T8ti2T] . Remarkably, 
in |2T] a good mathematical understanding of nested UDD (NUDD) is presented to 
suppress decoherence in arbitrary N-qubit system. NUDD scheme is based on mutually 
orthogonal operation set (MOOS). In addition, there exist more general DD schemes 
to protect a multi-qubit system, i.e. the periodic DD (PDD) and the symmetric DD 
(SDD) [?]. Both PDD and SDD are based on an appropriate DD group, and the two 
schems can be used to eliminate errors to the first and the second order in the Magnus 
expansion, respectively. 

In this paper, we aim to provide some guidelines in the choice of sequences to be 
applied in experiments on the suppression of decoherence for a multi-qubit system. To 
realize the goal, firstly we employ the concept of DD pulse-sequence construction as a 
filter-design problem [22] and derive the filter functions according to NUDD and SDD 
in the pure-dephasing spin-boson model with N qubits. The qubits embedded in 
two limiting cases, a common environment and separated environments, are considered. 
Furthermore, we discuss in detail the performances of NUDD and SDD for two qubits. 
In the case of a single qubit, it is shown that the UDD sequence outperforms the CPMG 
sequence for pure dephasing noises with sharp high-frequency cutoff while it performs 
slightly worse for soft cutoffs [231427] . As for a two-qubit system, it is also significant to 
discuss the performances of NUDD and SDD and see which performs better under the 
soft-cutoff and hard-cutoff baths. How to simulate the soft-cutoff and hard-cutoff baths? 
In terms of an experimental point, the noise used to examine the DD sequences, such 
as classical noise and ohmic spectrum, is in general governed by a power law spectrum 
l/w" (hereafter, a > 0) [28]. Varying the exponent of the power law of the spectrum 
allows us to simulate both baths with soft and hard cutoffs. On the other hand, how 
to test the performance of DD sequence? Exactly, a recent study provides us a novel 
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way, which is determined by the sequence itself only, to measure the performance on 
the suppression of decoherence when the system is subject to the noise scaling such as 
oc l/cj" [28j. The result shows that, except for the case of some small total number of 
pulses, SDD performs better than NUDD for the bath with the l/oo (soft-cutoff) noise 
and NUDD almost coincide with SDD but is still outperformed by SDD for the bath 
with the l/cj^ (hard-cutoff) noise. In addition, SDD helps to save the decoherence- 
free subspace while NUDD damages it. As the total pulse number increases, all filter 
functions for the same DD sequence have the same performance gradually in the two 
limiting cases mentioned above. 

Then a new analytical framework of filter-design perspective is employed to analyze 
some of the results observed [22]. The filter functions for NUDD and SDD in 
the coherence time regime and high-fidelity regime are compared, respectively. The 
definition of the two regimes and the differences between them are given in |22]. Both 
high-fidelity and coherence time regime are ranges of frequency for a certain filter 
function. Filter functions are designed to increase the suppression of errors in high- 
fidelity regime at short times, while the error probabilities of a few tens of percent are 
accepted in coherence time regime. From the comparison, we find that the function 
in the coherence time regime for SDD with the same total pulse number as NUDD 
peaks towards higher frequencies than that of NUDD, which is similar to the difference 
between CPMG and UDD. The functions in this regime can be used to demonstrate 
why SDD performs better than NUDD. In addition, the filter function for SDD in the 
high-fidelity regime has the same roffoff for various values of total pulse number while 
NUDD has different roffoffs for different values of total pulse number. The control pulses 
in the studies now are assumed to be ideal pulses, which is unphysical. Therefore, the 
control pulses with finite-width pulses are also considered in this paper. We find that 
the finite width affects the coherence time regime only for SDD sequence while it has 
influence on both coherence time regime and high-fidelity regime for NUDD sequences. 

The paper is organized as follows. In section 2, the dynamical decay of N qubits 
under dephasing is analyzed and its evolution operator with general DD strategy, SDD 
or NUDD, is derived in terms of a filter. In section 3, a two-qubit system, as an 
example, is disscussed in detail and the result about it is demonstrated by the filter- 
design perspective. The conclusion will be given in section 4. 

2. Decoherence suppression in a multi-qubit system 

2.1. Dynamical decay under dephasing 

We consider N qubits which do not interact directly among each other, while each qubit 
may have a different coupling to the bath modes p9|l30]. The total system Hamiltonian, 
in units of h, reads 



H = Hs + Hb + HsB, 



(1) 
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with 

N-l (j) 

Hs = Y.\-'^'^ (2) 

j=0 

-f^B = ^ ^kh\hk , (3) 

k 

N-l 

j=0 k 

where the first and second contribution Hs and Hb describe, respectively, the free 
evolution of the qubits and the environment, the third term Hsb describes a bihnear 
interaction between the two. Ua'' and a'^P are the transition frequency and the inversion 
operators of the jth qubit, h\ and hk are creation and annihilation operators for the kth 
field mode, which are characterized by a set of parameters {q"^\(jJk]- This information 
is encoded in the spectral density 

k 

The Hamiltonian Hsb in Eq-dl]), in interaction picture with respect to the free 
dynamics {Hs + Hb), is given as 

Af-l 

#5B(t) = $^ap)5(^)(t), (6) 

j=0 

with 

B^'\t) = J2 {ik^bU''-'' + ilk^rhe'''''')- (7) 

k 

The evolution is determined by the time-ordered unitary operator 

Uf{t)=fexp\ -I [ dsHsBis)}. (8) 







Under two standard assumptions: (i) the qubits and the environment are initially 
uncorrelated and (ii) the environment is initially in thermal equilibrium at temperature 
Te, the reduced density matrix of the multi-qubit system, in the standard product basis 
1 = {|0), . . . , - 1)} ( an N-digit binary notation, e.g., |6) = IIIO2) = | tU) ), 
can be written as 

psit) = TTB[Ufit)psiO)PBm}it)y (9) 
The element of this reduced density matrix is 



{m\ps{t)\n) = pmnif) 



Pmn(0)((K[„W)V-„(t)), (10) 
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with 

y±Jt)=fexp{±W dsY,{m,-n,)B^^\s)], (11) 

where irij G {0, 1} indicates the jth digit in the binary representation of the number 
m, for example, the 0th, 1st and 2nd digit in the binary representation of number 6, 
i.e.|6) = IIIO2), are 60 = 0,6i = 1,62 = 1, respectively. 
Thus we obtain 

Pmn(^) Pmn(0)Cm7i(^)? (12) 

with 

Cmnit) = (f expj / ds 5^(m, -n,)5(^)(s)}). (13) 

Note that if mj = uj the coherence function ( |T3l) is zero and then the diagonal 
elements are constant, as expected in a pure dephasing model. 

So far, the evolution without pulses has been analyzed. In the following subsection, 
we incorporate sequence of pulses into the dynamical evolution so that the decoherence 
of the system can be suppressed. 



2.2. Filter function formalism 

To suppress the decoherence in the pure-dephasing spin-boson model, we can apply a 
certain number Dj of ideal (5-shaped) vr pulses about the x axis on the jth qubit at 
times T^^'\ for d = 1, . . . , Dj, during the time interval from Tq'''' = to T^^^^-^ = T. Note 
that the pulses for different qubits may be applied at the same or different times (i.e., 
the pulse operators may be single-qubit or multi-qubit operations). After the pulses are 
incorporated into the free dynamical evolution, the Eq. fll3l) will be modified into [31j 



with 



Cmn(T) = (WmniT)), (14) 



T N-1 



W^n{T) = f exp I - 22 ^ ds [Y^i^j - n,)B'^^\s)F'^^\s)'\ } 



(15) 



D,-l 

F^'\t) = E (-l)'^(^ - ^i'V(TiJ, - t), (16) 

where the step function 6{t) is equal to 1 if t > and if t < 0. Here we discuss two 
limiting cases [29]: (i) Qubits feel a common reservoir, i.e. q^^^ = qk ■ In this extreme, 
the separations of the qubits are small compared to the wave length of the field modes, 
(ii) Qubits see independent reservoirs, i.e. q^^'' = qky Contrast to the former, in this 
case the qubits are so far apart from each other that each field mode couples only to 
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and 



(17) 



a single qubit. We use the P-representation for the thermal density matrix 
obtain the expectation value of Wmn 

Cmn{T) = exp I - ^ FmniujT)^ 

where the equal sign with a dot above ( i.e. =) represents c-number phase factors have 
been omitted, and the noise spectrum S{oj) is related to the spectral density J{ou) in 
Eq.® by 



S{uj) = J{uj) coth 
and the filter function [3ll[35] for q^f^^ = q^ reads 



(18) 



(cT) 



N-l 



j=0 



(19) 



while for q^^^ = qk^, we arrive at 



N-l 



T) = Y,\{^^-n,)\\f^\ujT)\ 

j=0 



(20) 



with the sampling function 



f^^\ioT) = -iuj 



i-T 

/ e-''''F^^\t)dt = l+(-l)^^+ie-'"^+2 V(-l)'^e- 
•^0 d=l 



-7^r,(2l) 



which encapsulates actually all information about an arbitrary sequence. The 
comparison between filter function f|T9|) and fl20|) clearly shows that the sum in Eq.f jTOj) 
can be understood as the interference of the sampling functions applied on different 
qubits, while the sum in Eq.( l20|) doesn't indicate this interference effect. A question is 
what different performances DD schemes have in these two limiting cases. This is one of 
the subjects of this paper. On the other hand, to measure the performance of different 
DD schemes, an approach proposed by Pasini et al [28]. can be used. The approach is 
suitable for power spectrum S{u) = ^(a > 0). We can write down the construction as 

So ^22) 



a+2 ' 



Then the decoherence function in Eq.( |T7I) gives 



[ du^^ Fmn{(^T) = So [ duj^—Fmn{(^T). 

Jo ^ Jo ^ 



Xmn 

(T) 

Substituting z = coT in the Eq.(!23l). we obtain 



with 



ya+2 



dz. 



(23) 

(24) 
(25) 
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The Eq.( 121|) reveals that the time dependence of XmniT) is a simple power of T. 
Meanwhile, the factor / is a quantity which expresses how well the pulse sequence 
protects a system against unwanted evolution. Note that we can simulate both baths 
with soft and hard cutoffs by varying the exponent of the power law of the spectrum. 
The smaller exponent a is, the softer the UV behavior of the power law spectrum is. 
Vice versa, the larger exponent a is, the harder the UV behavior of the power law 
spectrum is. The famous 1// noise corresponds to a = 1 in the notation. The case 
S{uj) (X Xjbj^ is experimentally relevant for ions in a Penning trap [25l|26]. This case 
corresponds to a = 4 in the above notation. The second question we want to report is 
which DD schemes perform better for softer and harder baths. 

We next describe two kinds of DD schemes, SDD and NUDD, for the model we 
consider and give their specific expressions of Eq.l pTi) . Then non-ideal control pulse 
with finite-width is also considered. 



2.2.1. SDD sequence SDD schemes is a group-based DD [7]. The element in the group 

{gi]i=o with order \G\, where g^ = Ig, is actually a base of linear unitary operators 
and the Hamiltonian, such as H in Eq.([l]), can be expressed as a linear combination of 
these bases. The building block of a group-based DD is realized by applying the pulses 
Pi = gi+igj separated by the same time delay At. Thus, the time evolution after a 
control cycle T = IGj At can be written as 

um = {gU-i^-'^'^'m-?! (^?fci-2e-^^^*^ici-2) ■ ■ ■ {gle-^'^'-'go) 

= fexp\-il Hg{s)ds], (26) 



where Hg{s) = g\e ''^^^gi, for s G {iAt, {i + l)At]. In the standard time-dependent 
perturbation theory formalism, the propagator is expanded up to the second order as 

um = i + J2h^ + H ^^^^ + ^ E + o{T'), (27) 

i i>j i 

where hi = —iAtHg{s). We shall say that fcth-order decoupling is achieved if the first 
k orders of the expansions of the propagator commutes with an arbitrary element of 
the group {gi}iJQ , and the first non-commuting term arises from the {k + l)th-order 
term. It was shown that |36], U{T) in Eq. fl27|) realizes the first-order decoupling, since 



we have the commuting correlation 



■splGI-l 7 



0, with the limiting T — )■ 



for i = 0, 1, . . . , |G| — 1, while to the second-order term we haven't this commuting 
correlation. This group-based cyclic sequence is referred to as PDD. Furthermore, we 
can obtain second-order decoupling via so-called SDD. The cycle becomes twice as long 
as PDD, T^^^ = 2T, and the time evolution operator is given by 

USDD^2T) = UU =1 + 2 J2 hi + ^ \ 2J2 hi] + 0{T^) (28) 
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where U = e'*°e^^ . . . e'^i'^i-^e'^i'^i-^ is mirror symmetric with U. In this case first-order 
and all even-order terms commute with gi, which generally makes SDD better than 
PDD. As for an N-qubit system [7], a possible choice of the group is G = cr^"'^}®^, 

P = x,y,z, j = 1,...,N, \G\ = 4^. However, if the relevant system-bath coupling 
is known to be linear in single-qubit operators, Hsb = "^^^^ ® ^ 

simplified group, G = {Is,^fSQ^o-^p^}, P = x,y,z, j = 1, . . . , N — 1, \G\ = 4. In 
particular, a purely decoherence coupling, like Hsb in Eq.(j4]), the decoupling group can 
be further simplified into G = {Is, ®f=oCx^}i \G\ = 2. Therefore, the SDD scheme 
for the system considered in this study is just to repeatedly apply the following pulse 
sequence: (fXfX)(XfXf)=fXffXf, where f and X denote a "pulse-free" period of a fixed 
duration and the matrix ^f=Qax \ respectively. Then the sampling function ( 12T]) can 
be specified as 

where Dj stands for the number of pulses applied on the jth qubit and is an even 
number, the same value for every qubit, that is, Dq = Di = ■ ■ ■ = -Dat-i. 

2.2.2. NUDD sequence NUDD schemes is a MOOS-based DD [2T]. A MOOS is 
defined as a set of operators which are unitary and Hermitian and have the property 
that each pair of elements either commutes or anticommutes. We can eliminate 
the effects of unwanted interactions between a system and its environment by the 
protection of this MOOS. A choice of the MOOS for our pure dephasing model is 

Now we describe the construction of NUDD for protecting {crx^}^S(^^. The 
construction consists of levels of control, the zeroth level, . . . , the (A^ — 2)th level and 
the (A^ — l)th level. First, in the (A^ — l)th level (the outermost level), Li\f_i operators 
of ai^ are applied at UDD timing 

Ti^_, = Tsin^ 2lXIi+2 ' In-i = ^,---,Ln-i , (30) 

between Tq = and T2,^_^-|_i = T. L^^i could be either odd or even. Then the free 
evolution in each interval is substituted by Ln^2 operators of ai^ applied at 

in each interval between Tij^_^fi = Ti^_^ and Ti^_^,Ln _2+i = ^U-i+i with Ln-2 being 
an even number. This is the (A^ — 2)th level. So on and so forth, the nth level 
of control is constructed by applying L„ operators of ai"^ in each interval between 

Tlf4-\,...,ln + l,ln ^Zjv_l,...,^n + 1 ~^ ( JV- 1 , ■ ■ ■ ,'n + 2 + l + 1 "^'jV- 1 1 • • -I'n + l ) 2L„+2 ' 

for ln = l, 
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with Ln being an even number. Thus, the samphng function (121 p for the (A^ — l)th, 
(A^ — 2)th and nth level can be specified as 

C-Dli^T) = e'^ Y: (-l)'-e^-^(^) (33) 

ijv-i=— Ljv-i — 1 



fZfo('^T) = E E (-l)"-"e^"-'^'. (34) 



-L 



N- 



Ljv-1 L„+i L„ 

„('t7~1 , „, \ ^ V— \ 'jV-1 >-'->'n + l V ^ , 'iV-l-----'ri + l inTT ^ 

At^^(u;r)= Yl ■■■ E ^ ' E (-1)^ ^ (35) 

/jv-i=0 i„+i=0 Z„=— L„ — 1 

where the evolution intervals r;^_^ = T,^_^+i - T;^_^ and rij^_j,...,i„+, = Ti^_^,...,i„+i+i - 
^/iv„i,...,/„+i- It should be noted that the nth nesting level of NUDD can be protected 
to the L„th order, while the overall protection of the MOOS is up to the Lth order 
with L = mm{Liy^i, Ln-2, ■ ■ ■ , Lq}. Athough NUDD scheme can protect a multi-qubit 
system in pure dephasing bath to a higher decoupling order than SDD, it remains a 
question whether NUDD offers improvement compared to SDD using the same total 
number of pulses and the same total pulse sequence duration T. We next discuss this 
question in detail. The reason for fixing the total time is that when T is set to a certain 
value, SDD has the same pulse interval for every qubit while NUDD has different pulse 
intervals for different qubits, i.e. intervals for NUDD may be larger or smaller than that 
for SDD. Therefore, evaluating which sequence it might favor is not easy if T is a fixed 
value, which ensures a relative fair comparison. 

2.2.3. Realistic pulses The above theoretical studies are under the assumption that the 

duration pulse time r^r = 0. When the pulses are much shorter than the other timescales 
of the system and the bath, it is a good approximation to treat them as infinitely short. 
But this assumption is unphysical since not all real system will satisfy the precondition. 
Here we consider a more realistic 7^ 0, and suppose that the interaction between 

system and bath is negligible during the application of each pulse. Then the general 
sampling function fl2T]) incorporated a nonzero r^r is modified into [25|[26] 

/(^■)'^(a;T) = 1 + {-if^+^e-'^T + 2 ^(-l)^e-^'^^^''"' cos(^), (36) 

d=i 

where Tjf^'"^ is the center time of the dth vr pulse. Therefore, the sampling fuctions (|29|) 

for SDD sequences can be modified into 

& = cos(^)/S., + [1 - cos(^)](l - e-^), (37) 
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and the sampling functions (|55]) ~ fl35p for NUDD can be rewritten as 

= cos(^^)/f^-^i + [1 - cos(^^)](l - e-n, (38) 

Ln-1 

C-nt = cosi^) fZ-n% + [1 - cos(^)] (1 - e"---.), (39) 

ijv-i=0 

fPuDD = ^os(^)fP^^^+[l-cos(^)] ... 5^ (1-e— '-.^■■■^W),(40) 

/]v-i=o ;„+i=o 

where ti^_^ and r/^_^,...,i„_^^, as before, are defined as ti^_^ = Ti^_^+i - Ti^_-^ and 
Tzjv_i,...,z„+i = ^/Ar_i,...,«„+i+i - ^Ziv-i,.-,«n+i- In this study, we also analyze whether the 
roles of for SDD and NUDD are the same or different. 



3. Numerical results and discussion 



For simplicity, we next restrict our attention to two-qubit system initially prepared in 
an arbitrary X two-qubit state, with only diagonal and anti-diagonal elements, under 
the stand product basis B = {|1) = | tt), |2) = | U), |3) = | it), |4) = | U)}- The 
MOOS for NUDD sequence is {ai^^}]=o = Wx\(t^x^}. From Eq.(l33])~(l35l), we can write 
down the sampling functions of this situation as 

e^ (-l)'ie^^-(^) (41) 

h=-Li-l 



fNUDoi^'^) - 



Cunn(-T)^±e^ f, (-l)'.e^-'#i', (42) 

li=0 lo=-Lo-l 

with 

Ti, = T,,+i - Ti, (43) 

where T^^+i and Tj^ can be obtained from Eq.( l30i) . and Li is an even number here. 
Accordingly, when the non-ideal pulses with finite width r^r are considered, the sampling 
functions can be modified into 

cos(^)](l-e-^), (44) 
cos(^)]|^(l-e— '0- (45) 



J NUDD ~ o ) J NUDD ' 



J NUDD — ^^°K 2 ) J NUDD ^ 
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On the other hand, the decouphng group for SDD sequence is G = {/s, ai'^Vi^''}, \G\ = 2. 
The samphng functions and the modified samphng functions read 

4ze'^ sin(^)sin2(|jf ) 

& = cos(i^^)/So + |1 - coB(i^)l(l - e-n, (47) 

with j = 0,1. Note that the pulse numbers Dq and Di are even numbers and equivalent. 
From Eq. ( IT9|) and ( 120|) , we have the filter functions for the case of a common reservoir 

F^4(u;T) = |/W(a;T) + /«(a;T)f, (48) 

F^,{coT) = \f^^\coT)-fW(coT)\\ (49) 
and the filter functions for the case of two separated reservoirs 

FUcoT) = FUcoT) = \f^'\coT)\'+\f^'\coT)\\ (50) 



Now, the performances of NUDD and SDD sequences in the two limiting cases, 
a common bath and two separated baths, are discussed. To ensure a relative fair 
comparison, in this study we assume both the total number of pulses and the total pulse 
sequence durations are the same for these two schemes, i.e. Lq + (Lq + l)Li = Dq + Di 
and T^^^^ = T^^^. In fi gure dl the factor / [2H] , see Eq. f l25p , is used to measure the 
performances of NUDD and SDD sequences for different filter functions F^^, F^^, F^^ 
and under the spectra S^u) oc 1/u and S{u!) oc The factors for different filter 

functions with the same sequence have obviously different performances if the number 
of pulses is small, while they almost coincide if the number of pulses becomes larger. In 
addition, the evaluation of the factor shows that SDD performs better than NUDD for 
the same filter function under S{uj) cc l/u, while under S{u) cc NUDD performs 
similarly to SDD but is still outperformed by SDD over the range of larger number of 
pulses. Note that there is a special dot {Dq,Di) = (4,4) for the SDD sequence with 
filter function F^^, where it performs abnormally and worse than NUDD for the same 
filter function. We will come back to this point below. It is also interesting to notice 
that the data for SDD with is always zero but it is not for NUDD sequence with 
F22. This means if the qubits are coupled to a common reservoir, SDD helps to protect 
the decoherence-free subspace while NUDD destroys it. 

Then, we employ the newly proposed approach [22] .including the figure layout, 
to discuss why SDD performs better than NUDD and this question can be clarified 
in the coherence time regime. The coherence time regime is a range of frequency 
for a filter function, in which the error probabilities of several tens of percent are 
accepted. The accumulated error leading to decoherence in this regime is very large. 
For clarity, we focus on the filter functions, F^^ and Fl^, and wrap the l/u"^ term into 
these functions, yielding modified filter functions F^^/u"^ and F^^/u"^. As shown in 
Eq.( l23|) . the decoherence function x(T) consists of the product of the spectrum S{u) 
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and the modified filter function. The decoherence function is minimum if the overlap 
between S{(jj) and the modified filter function is minimum. The power law spectrum 
S(u}) considered here is a strong low-frequency noise and the intensity decreases as 
frequency increases. Smaller values of modified filter functions in the low-frequency 
spectral regime will lead to smaller values of decoherence function x{T)i and hence 
smaller factor /, supposing T^^^^ = T^^^. In figure [2], the filter functions as a function 
of dimensionless angular frequency is shown for NUDD and SDD sequences with the 
pulse numbers considered in figure [H As expected, both NUDD and SDD serve as 
high-pass filters and increasing total pulse number corresponds to a shift in spectral 
peak towards higher frequencies. Thus figure [1] reveals that the factor I decreases with 
pulse number increasing. Furthermore, comparing figure [2l^a) with figure [2t^b) and figure 
M^c) with figure [2]^d), we observe the modified filter function for SDD peaks at higher 
frequencies than that for NUDD with the same total pulse number when this number 
becomes larger. Therefore, SDD performs better than NUDD when the pulse number is 
large, as shown in figure [H However, compared with SDD, the modified filter function 
F^^/u'^ for NUDD peaks at almost the same frequencies but has a smaller value when 
the pulse number is small, such as {Lo,Li) = (2,2) or {Do,Di) = (4,4). So it can be 
expected that NUDD performs better than SDD for filter function F^^ with this small 
pulse number, as shown in figure [H 

Here, the comparison of the filter functions between NUDD and SDD in another 
regime of interest, the high-fidelity regime, is also given. The high-fidelity regime, like 
the coherence time regime, is a range of frequency for a filter function, but in this regime 
the accumulated error bringing about decoherence is small. It turns up at short times 
t <^ T2, where T2 is the 1/e coherence time of a system. As before, we restrict our 
attention to the filter functions F^^ and F^^. The filter functions for the high-fidelity 
regime are shown graphically in figure [3l where the filter functions of NUDD and SDD 
are presented on a log-log plot for various values of (Lq, Li) or {Dq, Di). From figure[3]^a) 
and[3]^c), we find that NUDD sequence provides a filter function whose low-frequency 
rolloff increases from ^ 18dB/octave to ^ lOOdB/octave as pulse number increases 
from (Lo,I/i) = (2,2) to (Lo,I/i) = (16,16). By contrast, the low-frequency rolloff 
for SDD is approximately constant, ^ 18dB/octave, with each {Dq,Di). Here, rolloff 
is a term describing the steepness between the passband and the stopband of a filter 
function. The filter filters the noise in the stopband more efficiently, when the rolloff 
is steeper. Increasing rolloff means increasing the order of a filter and bringing the 
filter closer to the ideal response. If the rolloff is GrjdB /octave, the order of the filter is 
considered as rj. Then, as the total pulse number increases gradually, the filter order for 
NUDD becomes higher while the filter order for SDD stays about the same. Therefore, 
NUDD gives more fiexibility and performs better than SDD, in terms of suppressing 
the low-frequency noise. Note that the performances of NUDD and SDD are analogous 
to the performances of UDD and CPMG in the high-fidelity regime for a single qubit, 
respectively [22]. The coherence time regime is concerned with the decoherence function 
while the high-fidelity regime provides small contributions to decoherence, but the high- 
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fidelity regime draws more stringent scrutiny than the coherence time regime, because 
the high-fidehty regime is connected with predicted fault-tolerance error thresholds of 
Pth = 0.01% derived from quantum error correction and the maximum allowable error 
must not surpass pth in quantum computing applications [22]. The illustration about 
these two regimes here helps to understand what changes when the pulses have nonzero 
in the following. 

Now let us compare the finite-width filter function with the ideal filter function 
[22l[25] to find out what is affected by the finite-width pulses. In figure HJ we plot the 
ratio of the finite- width filter function F^^ (or Fl']^) to the ideal filter function F^^{ot F^^) 
for NUDD and SDD with a given total pulse number and here we set (Lq, Li) = (6, 6) 
or {Dq,Di) = (24,24). In these figures different traces correspond to different values 
of r^r in units of the total pulse duration time T. From figure ID^a) and ID^c), it can 
be seen that the filter function for NUDD is affected in both coherence time regime 
and high-fidelity regime, still there is a mid frequency range which is nearly unaffected. 
Moreover, as the pulse duration r^r increases, the influence on the filter function in the 
high-fidelity regime is more remarkable, and the mid frequency range becomes narrower. 
In comparison, figure Hl^b) and ll](c) show that the filter function for SDD is nearly 
unaffected in the high-fidelity regime while the magnitude of the changes in coherence 
time regime is larger than that for NUDD. Similarly, for SDD the range which is nearly 
unaffected becomes narrower with the increasing of pulse duration r,r- 

Fig.4 is ploted numerically, in which the sampling step is 0.0001 from 0.001 to 10 
and 100 from 10 to 10^. In the following, more details will be present. We first give out 
the analytical expression about the ratio of the finite-width filter function to the ideal 
titer function. For convenience, we focus on the case of two independent reservoirs and 
the expression for SDD reads 



j5 ^14:,SDD 

^SDD — -TTi 

^IA,SDD 



= cos (— ) + o.^TN cos(^) [ ■2(u.T . + 2 cos(— ) cos(a;T)] (51) 

It can be seen from Eq. flSTj) that there are singular points uT = AkDjir, with k 
being integral number and sin^(^) 7^ 0. On the other hand, the expression for NUDD 
is written as 



R 



NUDD 



^IA,NUDD 
^li,NUDD 



yjRe ) "I" yJim ) "I" yjRe ) + \Jlm ) /f-r,\ 

(/rp+(0=+(/S)^+(/a'^ ^ 
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where and f'^^}'^ ■, fj^'^ {n = 0, 1) are real parts and image parts of the ideal 

filter function and of the imperfect filter function, respectively. They are given by 

An) _ An) -An) /p-oN 
J NUDD — J Re ~^ ''J Im W'-'J 



An),r _ rynj,r ■ f(n),r /p. .x 

J NUDD — J Re ''J Im W^J 



with 



iRe = E (-l)'Hcos(— )cos[— cos(-^i— )] + sm(— )sm[— cos(-^^)]}(55) 



fil = E (-l)'^{-cos(— )sm[— cos(-^i— )]+sm(— )cos[— cos(-^i— )]}(56) 

h=-Li-l 

/jee =2^ 2^ (-1) {cos(^) cos[^ cos(-^-^)]+sm(^) sm[^ 



Zi=Oio=-io-l 
Li Lo 



fim = l^ 2^ (-1) {-cos(— -)sm[— -cos(— — -)]+sm(— -)cos[— -cos(-— — 

;i=oio=-i()-i ° ° 

/r = »s(^)/Ji' + 4sm^(i^)sm^(if ) (59) 
= ■:os(^^)/i;i + 2 8m^(^) 8m(^r) (60) 

Frc = cos(— +4sm (— ) 2^sm (— -) (61) 



2 ' 4 



Zi=0 



/i^" = cos(^)/iS + 2sin^(^) Esin(a;r,) (62) 

/i=0 

It is a tough task to list the singular points of Eq. (152|) one by one, but we can study 
this problem numerically. Figure 5 is the enlarged view of figure 4(c) and (d) within 
three ranges, i.e. from 250 to 1000, from 9.99 x 10^ to 1 x 10^ and from 3.999 x 10^ 
to 4 X 10^. From figure 5(a), we can see that the singular points appear periodically 
for the SDD sequency, as predicted by Eq. llSTI) . However, figure 5(b) indicates that no 
periodical singular points are obtained within the ranges considered. Therefore, the 
filter function for SDD is nearly unaffected in the high-fidelity regime, that is, before 
the first singular point (wT = ^D^n = OGn ^ 300, Dj = 24) appears. In addition, SDD 
is affected by imperfect control pulses in the coherence time regime, that is, after the 
first singular point {uT ^ 300) presents. The numerical results also suggest that in the 
coherence time regime the realistic filter function of NUDD is nearly equivalent to the 
ideal filter function. 



Filter-design perspective applying to dynamical decoupling of a multi-qubit system 15 
4. Conclusions 

In this paper we aim to provide some guidelines in the choice of sequences to be apphed in 
experiments on the suppression of decoherence for a multi-qubit system. We derived the 
filter functions for NUDD and SDD sequences with zero-width or finite-width pulses in 
the pure-dephasing spin-boson model. The filter design perspective was used to analyze 
the performances of these sequences. The performances of two qubits initially in X state 
were discussed in detail. 

The analysis shows that SDD outperforms NUDD for the bath with a soft cutoff 
while NUDD approaches SDD as the cutoff becomes harder. Second, if the qubits are 
coupled to a common reservoir, SDD helps to protect the dccoherence-free subspace 
while NUDD destroys it. Third, when the imperfect control pulses with finite width 
are considered, NUDD is affected in both the high-fidehty regime and coherence time 
regime while SDD is obviously affected in the coherence time regime only. Note that 
SDD scheme protects multiqubit systems using multi-qubit operations while NUDD do 
it using only single-qubit operations. Therefore, future work for SDD will analyze the 
effect of the precision, with which pulse location of each qubit in a multi-qubit operation 
is specified. In addition, the discussion for the case with more than two qubits will be 
given. 
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Figure 1. The factor / as a function of the number of pulses (Lo,Li) for NUDD 
sequence (squares) and {Dq.Di) for SDD sequence (circles). The pulse sequences are 
compared for different filter functions F{z) under the spectra (a) S{uj) oc l/w and (b) 
S{uj) oc 1/^4, see Eq.ipsl). 
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Figure 2. Modified filter functions, F^JiJ^ and FlJuP' ^ as a function of dimensionless 
angular frequency (wT) for NUDD and SDD sequences with the total pulse numbers 
studied in figure [TJ This figure is related to the coherencce time regime and its layout 
is referred to reference [55]. 
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Figure 3. Log-Log plot of the filter functions, F^/^ and Fl^, as a function of 
dimensionless angular frequency (wT) for NUDD and SDD sequences with the total 
pulse numbers studied in figure [TJ This figure is related to the high-fidelity regime and 
its layout is referred to reference P^ . 
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Figure 4. The ratio of the finite-width filter function F^l (or Fl'l) to the ideal 
filter function i^f4(or Fl^) for NUDD with (Lq, Li) = (6, 6) and SDD with (L>o, Di) = 
(24,24). In these figures different traces correspond to different values of Ttt in units 
of the total pulse duration time T . 
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Figure 5. Enlarged view of figure 4(c) and (d) within the dimensionless angular 
frequency ranges from 250 to 1000, from 9.99 x 10^ to 1 x 10^ and from 3.999 x 10^ to 
4 X 10^ for (a) SDD sequence and (b) NUDD sequence. The parameters are the same 
as that of figure 4. 



